Phonon-induced entanglement dynamics of two donor-based charge quantum bits 
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The entanglement dynamics of a pair of donor-based charge qubits is obtained in analytical form. The disen- 
tanglement is induced by off resonant scattering of acoustical phonons in the semiconductor host. According to 
our results a rather unusual recovery of entanglement occurs that depends on the geometrical configuration of 
the qubits. In addition, for large times a non-vanishing stationary entanglement is predicted. For the cases of one 
and two initial excitations a simple kinetic interpretation allows for an adequate analysis of the observed dynam- 
ics. Our results also reveal a direct relation between the disentanglement rate and the inter-donor decoherence 
rates. 

PACS numbers: 12.40.Nn, 11. 55. Jy, 05.20.-y, 05.70.Fh. 



I. INTRODUCTION 

Quantum information processing promises highly efficient 
solutions to cryptographic problems and exhaustive database 
search, that outperform the best known algorithms with clas- 
sical computers jlj. Such quantum algorithms rely on the ca- 
pability to process correlations among quantum subsystems. 
The quantum part of these correlations is denoted as entangle- 
ment. Most commonly, the subsystems are identified as two- 
level systems, representing a quantum counterpart of classical 
bits, called quantum bits (qubits). Among the well known ap- 
plications of quantum information are quantum teleportation 
yt], superdense coding yj], and secure distribution of crypto- 
graphic keys Jit. 

Entanglement, being the key ingredient of such applica- 
tions, at the same time is highly fragile and can be easily de- 
teriorated during state preparation, addressing and control of 
individual qubits, and final readout. Hence, it is necessary to 
encode quantum information in physical systems, where en- 
tanglement is protected or can be preserved in a robust way 
from ambient effects. Prominent examples of physical imple- 
mentations are trapped ions JH, nuclear magnetic resonance 
|0], atoms in cavities 0], quantum dotsjUt], semiconductor 
impurities jsjl, superconducting qubits IllOtl - and impurities in 
diamond ifTlll . 

Among the solid-state implementations Il2l - [l4ll impurities 
embedded in a semiconductor substrate offer the advantage of 
comparably easy scaling and production due to highly devel- 
oped fabrication techniques. The encoding of qubits in the 
charge degrees of freedom of pairs of donor sites allows for 
the addressing of individual qubits by metallic gates lfl5i [lcjll . 
Tunneling of the electron between donor sites together with 
Coulomb repulsion of electrons bound in neighboring qubits 
have been shown to allow for the realization of a CNOT gate 
lfl5i[l7i[T8ll . Such a gate is the elemental building block of any 
quantum algorithm. Typical coherence times of such systems 
are limited by phonon scattering to the order of lps Ill9l42lll . 
Therefore, the state preparation process must be faster and one 
may expect that entanglement is rapidly lost on this time scale. 

In this work we show that to a large extent the entangle- 
ment survives beyond this time scale and therefore suffers 



further degradation only from other sources of decoherence, 
i.e. charge fluctuations on the control electrodes. The disen- 
tanglement dynamics is obtained in analytical form showing 
non-Markovian features, similar to the decoherence dynamics 
of a single qubit Il20l422ll . Furthermore, it is shown that the 
disentanglement rate is directly related to inter-donor deco- 
herence rates for the cases of one and two initial excitations. 
The structure of this rate can be explained by a simple kinetic 
interpretation that allows for the determination of the disen- 
tanglement rate from the geometry of the constituent donor 
sites. 

The paper is organized as follows: In Sec. [IT]the dynamics 
of a general N qubit system subject to off-resonant scatter- 
ing of acoustical phonons is derived. Using these results the 
dynamics of entanglement between two qubits is analytically 
obtained in Sec. Ill for the cases of one and two excitations. 
Finally, in Sec. IV we present a summary and conclusions. 



n. PHONON-INDUCED DEPHASING DYNAMICS OF 
QUBITS 

In this chapter we deduce the dynamics of the reduced den- 
sity operator of the qubit system, which is induced by off- 
resonant scattering of acoustical phonons within the semicon- 
ductor material. It is assumed here that 2N donor sites are 
present that conform N qubits, see Fig. [T] The results will 
later be specialized to the case of N = 2 qubits. The dynam- 
ics of a single qubit conformed by two donor sites has been 
derived recently by us in Ref. Il20ll . The present work goes 
beyond the case of two donors, which implies a more com- 
plex dynamics together with the possibility of studying the 
entanglement between pairs of qubits. 

We consider the situation where a single electron is con- 
fined to each pair of donors, conforming the N qubits being 
centered at the positions r&, where b = l,...,N labels the 
qubit under consideration. In addition, we assume that the 
distances between these qubits are much larger than the inter- 
donor vectors d/, within the qubits, i.e. |rj, — ry \ 3> \&b"\- 
Under these circumstances the tunneling of electrons between 
different qubits can be safely neglected. Considering the in- 
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Figure 1: Geometrical setup of N qubits formed by 2N donor sites 
in a semiconductor material. Blue and red donor sites correspond to 
qubit states m b = — i and = +i, respectively, where b = 1 , . . . , N 
denotes the qubit. Each qubit b has an inter-site distance d;, pointing 
from site mj, = — \ to = + j, and a position vector r/, pointing 
towards the center of the qubit. 



teraction of the qubit charges with acoustical phonons in the 
semiconductor substrate at room temperature or below, off- 
resonant phonon scattering is the main source of decoherence 
in the electron dynamics [20]. This type of electron-phonon 
interaction is accurately described by the spin-boson Hamil- 
tonian 
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where the individual interaction rate of the m^th donor site 
(m/, = ±1/2) within qubit b reads 
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Here a k are the annihilation operators of longitudinal acousti- 
cal phonon with dispersion relation (Ot = sk, s being the sound 
speed. Moreover, Mq and D are the mass within a unit cell of 
the semiconductor and the deformation constant, respectively. 
Each donor site is described by a s-wave ground-state with 
Bohr radius a b m , and the electronic transition of each qubit b 



is 



generated by the pseudo spin-| operator S/,, where 



^=2 



d b 



(3) 

with \r b ztdb/2) being the states with the electron being lo- 
calized at the corresponding donor site. In what follows we 
assume that during the free evolution of the system, tunneling 
is inhibited by either an applied potential barrier between the 
donor sites, or due to a strong bias between the qubit levels, 
1 0>i| 3> | | , which can be provided for by the application of 
a DC electric field. 

Following the same steps as in Ref. lEoll . the Hamilto- 
nian Eq.(Q]) can be diagonalized, to obtain the complete set 
of eigenstates as displaced number states 

\ E {m b },{N k} ) = \{m b }) ®# ({« W} , k }) \{N k }), (4) 
with eigenenergies 



E {m b }.{N k } = hY, ( °l> m l> + Yj l0}kN ^ 
b k 

Here the multi-mode displacement operator reads 
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and |{A^k}} = Ilk l-^k) are multi-mode number states of the 
acoustic phonons and the displacement amplitude reads 
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where a b k = g b ,k/ ®k- Furthermore, the quantum state of the 
qubits is encoded in the register state 



\{m b }) = | mi) (8 |m 2 ) ® • • • ® \m N ), 



(8) 



where each of the qubits can be in states \m b = ±j) (b = 

1,..., AO- 
Given the eigenstates (01 and eigenenergies (O, the gen- 
eral solution of the reduced density operator of the qubits, i.e. 
traced over the phonons, reads 

Pit) = £ \{m b })({s b }\ £ E P{m b },{N x };{ Sb }M) 

{«*}.{«*} \Nk}{N k } 

x({Nl}\D({a {sb] . k })fr ({« W} , k }) \{N k }) 
it 



x exp 
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Here the matrix elements of the initial density operator of the 
complete electron-phonon system in the basis of the energy 
eigenstates are 

P{m b },{N k h{s„UK} = ^{'"*}>{A'k}lp( )l £ K}.{Afa^ (10) 

These density matrix elements are determined by the state- 
preparation process that is utilized to set up the initial entan- 
glement between the qubits. A generic state-preparation pro- 
cess can be described as follows: 
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We assume that initially the system is at low enough tem- 
perature, IcqT <C fia> b , such that it has relaxed completely into 
a state where all the qubits are in their lowest energy states, 
m b = — 1/2 (b = 1, . . . ,A0, and coexist in thermal equilibrium 
with the phonons in the substrate. Starting from this state 
the system undergoes a state preparation process in which the 
qubits can be coherently transferred into a superposition of 
ground and excited states, without affecting the quantum state 
of the phonons. The probability amplitudes to transfer from 
the initial qubit state { — i, — j,..., — ^} to the qubit states 
{m b }, shall be denoted by y/\ mb }- The corresponding tran- 
sition can be generated by the application of the operator 



S +,{**} - U b 



(11) 



where f b is the identity operator in the Hilbert space of qubit 
b. Thus, the matrix elements of the prepared initial state (TlOb 
become 



P{m b },{N k };{s b },{K} 



(12) 



where 

P T= E ^J^} [^{-1/2 -1/2,...},{A^}> <^{-l/2,-l/2,...},{JV k } I 
W 



is the initial thermal state with the phonon statistics 
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with Z satisfying £{jvjj ^{Ay = 1. and Pk = T«D k /(k B T). 

We note that the corresponding prepared initial reduced 
density operator of the qubits is pure, i.e. is of the form 
|i//'(0))(v(0)| with the qubit state being the sought superpo- 
sition 



|V(0)> = E ¥{m b }\{m})- 

i m b) 



(15) 



Thus, the presence of the phonons does not prevent a coherent 
preparation of the initial qubit state. However, the initial state 
preparation takes the electron-phonon system out of thermal 
equilibrium so that the decoherence of the quantum state of 
the qubits cannot be treated by the usual method of spectral 
functions l23ll . 

The generic state preparation, as described above, can be 
implemented for example by switching the gate voltages so 
that the on-site energies within a qubit cross each other in- 
ducing Landau-Zener transitions. At the end of the process, 
the system will be in a coherent superposition of all possible 



states of the system depending on how fast its two-level com- 
ponents where driven across the level crossing. Alternatively, 
it may be implemented by time controlled tunneling and em- 
ploying the Coulomb repulsion between neighboring qubits to 
generate entangled qubit states, as proposed for implementing 
CNOT gates (SoB- Furthermore, it could also be im- 
plemented by THz Raman transitions between the donor sites 

We note that the matrix elements of the initial reduced elec- 
tronic density operator are obtained from Eq. (0 as 

{{m b }\p s (0)\{s h }} = £ £ Pw,ft);( Sl ),«) 
WW 

with the non-diagonal elements containing the Franck- 
Condon type transition amplitudes 



({K}\D({a [sbhk })fr ({« W} , k }) \{N k }). 



(17) 



These factors are overlap integrals of two displaced phonon 
number states with displacements OC^ Sb y^ and tt^jjt, respec- 
tively. Their presence is due to the fact that the initial den- 
sity matrix of the complete electron-phonon system ( fTOb is 
in the basis of the energy eigenstates (@), whereas the ma- 
trix elements of Eq. ( [ToT l are in the basis of the product 
states that differ by the displacement of the 

phonons. 

Consistent with a dephasing model, the diagonal elements 
follow from Eq. (0 as invariants: 



({m}\ps(t)\{m b }) = ({m b }\ Ps (0)\{m b }). 



(18) 



However, this dephasing — being induced by off-resonant 
scattering of acoustical phonons — modifies the time evolu- 
tion of the off-diagonal density matrix elements as 



({m}\p s (t)\{s b }) 
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It can be observed in Eq. ( fT9l ) that apart from the free oscil- 
lation with angular frequency (Oq, a dephasing is induced by 
differing phonon numbers in combination with the presence 
of the Franck-Condon factor- 
To further evaluate the dephasing of the off-diagonal den- 
sity matrix elements of the qubits, we insert the initial com- 
plete density matrix elements ( fTOb together with Eq. ( fT3l into 
Eq. $1% . From this we obtain 
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X Xf {s b },{KY'{-1 -2}..{M k }A^b}4K}'{"y^'{ N k}-f{m b }.{N k };{-^,...-^},{M k }- 



(20) 



Employing the thermal phonon statistics ( fPfl i the sum over the phonon numbers in Eq. ( |20l can be rewritten as a trace, which 
leaves us with 
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where we defined the time-dependent phonon displacement 
amplitude 0!{ mA }, k (f) = 0£{ mij }. k exp(— iCO^t ). The displacement 
operators in Eq. (l2Tb can be combined to obtain 



xZ^Tr 



(22) 



with the time dependent phase being 



A {mfcH . s , } (0 = 3la W}lK}lk 5 { * mi}){ii})k - 1) (23) 



Here we have defined sum and difference displacements: 

«{m A M.v 6 },k = a{m i ,},k + «{ Si ,},k-2a { _i _i }jk , (24) 

5{m fc },{ Sfc },k = «{m fc },k-«{^},k- (25) 



The trace in Eq. ( 1221 represents a thermal average that can 
be evaluated in phase space to obtain 

({ mb }\p s (t)\M) = W{mb} r {Sb f-' u ^^ 

xexpj-^ dt'T {mbUsi} {t')\ , (26) 



where the decoherence rate of the qubit state is defined as 

T {m b },{s b }(t) = TV (')• (27) 

b,y 

Whereas this rate depends on the quantum numbers of the 
density matrix element under consideration, the relation of 
pairs of bits is governed by the inter-bit decorrelation rate, 
given by 

7b.b'(t) =4 Y J Yj mbS b'yVi a bM b , a b',s b nh,m b -,b',.s b ,) -(28) 

m b s b , 

This rate in turn depends via the inter-donor decoherence rate 
y(t;a,a\ I) on the Bohr radii of the four donor sites of the two 
qubits and on the six possible distances between these four 
donors, 



h,m b ;b\s y = \{*b+m b Ab)-{ry+syAy)\. 



(29) 



Moreover, the inter-donor decoherence rate is obtained from 
Eqs. d22l - d29|> as 



y(t;a,a,l) = Ft 
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In this expression the temperature dependent rate reads Tj = all donor sites and Nb is the number of unit cells within the 
(Qb(T /7b), where the convenient temperature scale is chosen average Bohr volume a^. 
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l,...,N), the inter-donor decoherence rate (130t approaches 
Here ©b = 2ns/a-& with «b being the average Bohr radius of the form y(t;a,a',l) — > y(t;l) with 
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Figure 2: Dimensionless and temperature independent inter-site de- 
coherence rate j(Iu)/Tt as a function of the traveled distance of a 
phonon in units of the Bohr radius, st/a%. 



We note that in the limit of vanishing distance between donor 
sites, / — > 0, this function becomes 



y(f;0) 



st 



St 



-2.vf/«B 



1 

K a B J 2 \a-B, 

(33) 

This is a function peaked at t ~ a% /s, i.e. at the time a phonon 
needs to travel the distance of one Bohr radius, see Fig. |2](red 
curve). Different from this special case (I —5- 0) for / > the 
inter-donor decoherence rate is peaked at the time t ~ //* that 
is required for a phonon to travel the distance Z, see Fig. [2] 
(green and blue curves). 



m. DYNAMICS OF ENTANGLEMENT BETWEEN TWO 
QUBITS 

In the following we will discuss the special case of two 
qubits being present in the semiconductor system, i.e. N = 2. 
In this case the state of the qubits with b =1,2 lives 

in a four-dimensional Hilbert space and the entanglement of 
this bipartite system can be described by the concurrence C 
ll25ll . This quantity is a measure of entanglement bounded 
in the range between zero and one, with the maximum en- 
tanglement corresponding to unit concurrence. On the other 
hand, for separable states, lacking any entanglement, the con- 
currence is zero. 

Furthermore, due to the pure dephasing effect of the phonon 
scattering, no transitions of the qubits are induced. Therefore, 
the phonon scattering will transform an initial general super- 
position state 



|V(0)> 



L V / {m l .m 2 }\{mi,m 2 }), (34) 
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into a non-pure statistical mixture of only those states that ini- 
tially already existed. Thus, during the time evolution the den- 



sity operator stays within the Hilbert sub space defined by the 
initial state. This feature allows us to separately treat the two 
prominent cases of having initially one or two "excitations", 
respectively. 



A. Case of one excitation 

Assume the initial state of the two qubits to be of the form 

Mo)) = Y { i,-^ } \{k-k}) + V{-^}\{-hh})- C 35 ) 

Since only one of the qubits is in its "excited" state, this 
superposition is usually denoted as the "one excitation" 
case. Choosing the basis vectors of the bipartite system as 

^}>J{| ,-|}>,|{4>i}>>l{-2 -2>>}> following Eq. 
(|2oT the time-dependent density matrix can be written as 



Ps(0 = 
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(36) 

For this particular form of the density matrix, the concurrence 
simplifies to 



C(t)=2 p.. 



(37) 



Inserting the corresponding density matrix element from Eq. 
into Eq. (|37| |, the concurrence results as 



C(0=2Vp(l-p)exp 



dt'T 



1 1 1 



K-h 



, (38) 



where p = \y^i _^j| 2 is the initial probability for the two 

qubits being in state ({5,— !})■ The maximum initial con- 
currence is of course obtained for equal weights, p = 1/2, of 
the two constituent states. 

From Eq. d38i l it becomes apparent that the decoherence 
rate of the state of the qubits acts as disentanglement rate. 
This rate is shown in Fig. [3]for qubits with d = IOab, an inter- 
qubit distance of 20a%, and a relative angle of 45°, see inset 
of Fig. [3] It shows a series of alternating maxima and minima 
at increasing times. The principal positive maximum at the 
beginning occurs at t « /s, which is the time needed by the 
phonon to travel within a donor site. This is the main source 
of disentanglement. The times of the subsequent extrema can 
be identified as the travel times between pairs of donor sites, 
as indicated in the inset of Fig. [3] Whereas the positive max- 
ima destroy, the negative minima restore the entanglement be- 
tween the qubits. The mapping of which phonon path between 
donor sites leads to positive or negative extrema in the disen- 
tanglement rate can be established as follows: 

The concurrence is given by the modulus of the den- 
sity matrix element 0,1 i w 1 i~,(t) that describes the 

12'~!/'1 I'll v ' 

time-dependent correlation between states |{ — 5,5}) and 
— A}). Correlations between these states can only be cre- 
ated when a phonon travels between a donor site occupied by 
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Figure 3: Dimensionless disentanglement rate Tri i^(t)/rx 

as a function of the dimensionless time st/a^ for the case of a sin- 
gle excitation. The inset shows the geometrical configuration of the 
donor sites: d\ = d.2 = IOqb. I r l — r l\ = 20ab, 45° angle between 
qubit axes. The rounded lengths of the inter-donor distances in units 
of a B are: l/a B = 10 (1), 16.53 (2), 18.6 (3), 23.6 (4), 25.0 (5). 



one state to another donor site that is occupied by the other 
state. The corresponding site occupations of each of the states 
involved are indicated in the inset of Fig. [3] by black and 
white colors, respectively. The sites of each qubit are marked 
by blue for = — A and red for ntj, = + A . With this color 
scheme, the creation of correlations mediated by sound waves 
is produced via phonon travels between a black and a white 
donor site. 

These phonon travels are: the passage within the individual 
qubits for the distance lOa^, that produces in Fig. [3]the nega- 
tive minimum 1, and the passages between blue and red sites 
at distances ps 16.5«b and ps 23.6«b, that produce the minima 
2 and 4. All the other phonon passages produce decorrela- 
tion and destroy the entanglement at distances ps 18.6«b an d 
« 25.0ab, which generate the positive maxima 3 and 5 in Fig. 

a 

We note that diminishing the 45° angle between the qubit 
axes results in lengths 2 and 4, and lengths 3 and 5, re- 
spectively, becoming progressively comparable. In the lim- 
iting case of two collinear qubits, i.e. 0° angle, these pairs 
of lengths are identical so that as a consequence the disen- 
tanglement rate shows only two positive and two negative 
peaks, as shown in Fig. |4] (a). On the other hand, as one 
approaches the limiting angle of 90°, i.e. the CNOT configu- 
ration [Tit], the extrema first sparse and finally peaks 
2 and 4 cancel peaks 3 and 5, respectively. As a result, the 
disentanglement rate shows only the principal positive peak at 
t p» gb /s, see Fig. H](b). 

The evolution of the concurrence for the case of 45° be- 
tween qubit axes is shown in Fig. [5] It can be seen that 
a stationary and non-vanishing value of the concurrence is 
reached for large times. Moreover, the temperature depen- 
dence indicates an only minor loss of entanglement at temper- 
atures T /7b < 0.01. Given that for P impurities embedded in 




Figure 4: Dimensionless disentanglement rate Tj-i _iir_i i}(0/Fr 

as a function of the dimensionless time st/a^ for the case of a single 
excitation. Same parameters as in Fig. 3, but for collinear (a) and 
perpendicular (b) qubits. 



a Si substrate the characteristic temperature is of the order of 
Tg ~ 300K, this case corresponds to liquid He temperatures. 



B. Case of two excitations 



The other prominent case is that of initial two excitations. 
This case is described by the initial state being of the form 



\Y(0)) = Y { ^J{h l 2}) + y { -l-^ } \{- 1 2 ( 39 ) 
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Figure 5: Concurrence as a function of the dimensionless time st/a B 
for the case of a single excitation for temperatures T/T B = 0.01 
(solid curve), 0.05 (dotted curve), 0.1 (dashed curve). Other param- 
eters same as in Fig. 3 with p = 1/2. 



In the same standard basis as before, the time-dependent den- 
sity matrix results then as 

/ P { - M - } (0 o o p { i,i m _i,_i } (0 \ 



Pv>- 

Vp { *- } , m ,_i } (0 o o P { _i_. M _._i } wj 

(40) 

Also for this case the concurrence simplifies to a simple ex- 
pression, given by 



P { ^},M,-^}( f ) 



C(t)=2 

which, after insertion of Eq. (), becomes 



(41) 



C(t)=2y/p(\-p)exp 



dt'Y, iiir i 



(0 



(42) 



where now p = \ y/^ i i ^ | 2 is the probability for the two qubits 
being "excited". From the change of indices it can be easily 
seen that, apart from the first positive and first negative peak, 
the signs of the peaks of the disentanglement rate are reversed 
as compared to the corresponding one excitation case, see Fig. 

m 



IV. SUMMARY AND OUTLOOK 

In this paper we considered the time evolution of entan- 
glement in donor-based charge qubits that is induced by off- 



resonant scattering with acoustical phonons. We showed 
that this system can be solved analytically and that a non- 
Markovian behavior emerges with negative disentanglement 
rates, leading to non-mono tonic disentanglement in time. 
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Figure 6: Dimensionless disentanglement rate F, i _i}(0AV 
as a function of the dimensionless time st/a B for the case of two 
excitations. The inset shows the geometrical configuration of the 
donor sites: d\ = di = \0a B , l r l — r l\ = 20«b> 45° angle between 
qubit axes. The rounded lengths of the inter-donor distances in units 
of a B are: l/a B = 10 (1), 16.53 (2), 18.6 (3), 23.6 (4), 25.0 (5). 



Moreover, for the cases of one and two initial excitations the 
disentanglement rate is proportional to the decoherence rate of 
the two-qubit state. In both cases the concurrence attains a sta- 
tionary and non-vanishing value at large times, which means 
that phonon scattering does not completely destroy the entan- 
glement of the initially prepared two-qubit state. 

The choice of the geometry of the donor sites determines 
the features of the concurrence as a time-dependent function. 
These features can be understood by a simple kinetic inter- 
pretation of phonon travels among the donor sites. In this 
work we focused on the cases of initially one and two excita- 
tions. However, we believe that also particular superpositions 
of both cases may be treated within this framework. Further- 
more, our model includes already the case of N > 2 qubits, 
where a trace over N — 2 qubits would be required to obtain 
the entanglement between a selected pair of qubits. This will 
be subject of future work. 
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